The Rotation Of The Deep Solar Layers 
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ABSTRACT 

From the analysis of low-order GOLF+MDI sectoral modes {£ < 3,6 < n < 15,|m| = £) and 
LOWL data {£ > 3), we derive the radial rotation profile assuming no latitudinal dependance 
in the solar core. These low-order acoustic modes contain the most statistically significant in- 
formation about rotation of the deepest solar layers and should be least influenced by internal 
variability associated with the solar dynamo. After correction of the sectoral splittings for their 
contamination by the rotation of the higher latitudes, we obtain a flat rotation profile down to 
0.2 Rq. 

Subject headings: Sun: lielioseismology, rotation — Instruments: GOLF, LOWL, MDI 



1. Introduction 

Helioseismologists use the oscillations of acous- 
tic waves that propagate inside the Sun to infer 
its rotation profile. Due to the solar rotation (and 
magnetic fields), the frequencies of two modes of 
the same degree £ and radial order n, but with 
different azimuthal orders m, are separated by a 
small amount referred to as splitting (Aujinm)- 

The understanding of the angular momentum 
redistribution in the deep interior requires a pre- 
cise derivation of the rotation profile below the 
convection zone, but this profile is still a matter of 
debate. In the past an increase of the rotation rate 
near the core has long been favored (e.g. Lazrek 
et al. 1996), but recent results by Chaplin et al. 
(1999) and Eff-Darwich, Korzennik, & Jimenez- 
Reyes (2002) favor a slight decrease. Corbard et 
al. (1997), using LOWL data, and Ehgamberdiev 
et al. (2001), using IRIS data, also derived such a 
decrease. 

Here we focus on the rotation in the solar core 
derived from the GOLF (Global Oscillations at 
Low Frequency) and MDI (Michelson Dopplcr Im- 



ager) instruments, supplemented for modes with 
£ > 3 hy LOWL data. These data sets and the 
interest of the low-order modes are presented in 
sections 2 and 3 respectively. In section 4, the 
method utilized to extract the rotational splittings 
is discussed. In section 5 we discuss the rotation 
profile. We conclude in section 6. 

2. The Data 

We use data from GOLF/SoHO and MDI/SoHO 
for the modes £ < 3. GOLF detects the global so- 
lar oscillations by the Doppler shift they produce 
at the surface on the sodium lines. It is specif- 
ically designed to be more sensitive to the low- 
degree modes, and only detects those with £ + m 
equal to an even number. We utilize the 2034-day 
long GOLF series (starting in April 1996) for the 
GOLF-alone splittings. MDI detects solar surface 
velocity from a nickel spectral line and resolves the 
solar surface. To produce integrated-disk tempo- 
ral series sensitive to the low-£ modes, we apply 
two masks (gaussian, and gaussian zero mean) 
to the 2243-day long MDI series (Henney et al. 
1999). These two series have been properly com- 
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bined with GOLF 2243-day long series to increase 
the signal-to-noise ratio of the \ow-£ modes, and 
to reduce the leakage of higher-^ modes present on 
the MDI spectrum (Garcia ct al. 2003). The use 
of the MDI-I-GOLF time series allows for the pre- 
cise measurement of lower-n modes. This point of 
interest is discussed further in section 3. Wc sup- 
plemented these GOLF-alone and MDI-hGOLF 
data with LOWL rotational splittings for £ > 3. 
We used six LOWL series of 1 year, from February 
1994 to February 2000. The splittings have been 
extracted for each year but the final result uses 
the 6 years statistics (Jimenez- Reyes 2001). 

3. Information Value of Low- Order Modes 

The derivation of the rotation profile in the ra- 
diative zone requires high-precision splittings. For 
instance, only ~ 3.5% of the splitting (406 nHz) 
of the mode £=ln = 9m=l(ed, 1472.85 /zHz) 
is due to the rotation below 0.2 Rq (Couvidat 
2002). Thus, we need to derive Awi^g^i with a 
precision better than 14.2 nHz to obtain informa- 
tion on the core. The error associated with this 
splitting is ~ ±5.3 nHz (see Table 1). It contains 
more information than the mode £ = 1 n = 20 
(at 2963.43 /j,Hz), for which 6.5% of its splitting 
(401 nHz) is due to the core, i.e. 26 nHz, but the 
uncertainty is 34.6 nHz. This point is illustrated 
by the left panel of Fig 1. Even though the in- 
ner turning point is closer to the core as n in- 
creases, the outer turning point is also closer to 
the surface, where stochastic excitation and solar 
cycle effects become more important. Therefore, 
the most favorable trade-off between sensitivity to 
the core rotation and uncertainty in the rotational 
splitting determination occurs for \ow-£ \ow-n p 
modes. Consequently we use only these modes in 
the present analysis. 

4. The Rotational Splittings 

The splittings are obtained by fitting the com- 
ponents of each mode to asymmetric Lorentzian 
profiles (Nigam & Kosovichev 1998) using a max- 
imum likelihood code based on Appourchaux, Gi- 
zon, & Rabello-Soares (1998). The modes £ = 
Sz 2 are fitted together, likewise the modes £ = 1 
& 3. For the GOLF-alone series, we use the peri- 
odogram of the data as the power spectrum, ex- 
cept for the modes £=ln = 7&8, and £ = 2 




Fig. 1. — Left panel: Ratio of the splitting due to 
the layers below 0.2 Rq, to the uncertainty on the 
total splitting for = 1 (solid line), £ ^ 2 (dashed 
line), and £ = 3 (dash-dotted line). Right panel: 
synodic splittings for n < 15 used for the inver- 
sion of the solar rotation, obtained with GOLF 
and MDI, as a function of the inner turning point 
for ^ = 1 (crosses), £ = 2 (stars), and £ = 3 (dia- 
monds). 

n = 8 & 9, for which we use a multitaper spec- 
tral estimator (e.g. Percival & Walden 1993). For 

the combined MDI-I-GOLF series we again use 
the multitaper. We restricted the search of the 
splittings to modes with a frequency less than 
~ 3000 /xHz, and use for the inversion only the 
modes with n < 15. This way we limit the im- 
pact of the solar cycle and surface effects. This 
approach is consistent with the results of Chap- 
lin et al. (2001) who demonstrated that the split- 
tings are strongly biased at higher frequencies, and 
Chaplin et al. (2003), who showed that the solar 
cycle does not affect the splittings of the low-n 
£ = 2 modes. 

The free parameters of the fitted profiles for 
each mode are: the full width at half maximum 
(FWHM, the same for all the components of a 
mode), the amplitude (the relative amplitudes of 
the different components are fixed empirically), 
the asymmetry (the same for each component), 
the central frequency, the noise level (assumed to 
be constant in the fit window of 40 /iHz), and the 
splittings. The latter parameters arc sectoral val- 
ues {\m\ = £). We restrict the fits to frequencies 
< 2000 - 2200 AiHz on the MDI+GOLF data, be- 
cause of the growing infiuence of higher-£ modes in 
the fit window (for instance, £ = f> modes perturb 
the ^ = 2 profiles). 

In Table 1, we list the synodic sectoral split- 
tings and their uncertainties. The rotational split- 
tings common to the MDI-^GOLF and GOLF- 
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alone datasets generally agree within the Icr un- 
certainty. Right panel of Fig 1 presents these split- 
tings as a function of the inner turning point. This 
plot emphasizes the high quality of the SoHO data 
after 5 years of observation in comparison with 
previous data sets. 

5. The Solar Rotation Profile 

The ID inversion of sectoral splittings pro- 
vides information about the rotation mainly along 
the solar equator, but with contamination by the 
higher latitudes. We applied a ID MOLA inver- 
sion technique (Multiplicative Optimally Local- 
ized Averages, see e.g. Corbard et al. 1998). The 
MDI-I-GOLF and GOLF-alone splittings (when 
the former are unavailable) were completed for 
i > 3 with splittings from LOWL. The direct in- 
version of the data leads to a rigid rotation from 
the base of the convective zone at 0.713 Rq down 
to ~ 0.35 ii© (see upper panel of Fig 2). Below 
there is a decrease in the rotation rate of the solar 
equator, as mentioned previously by several au- 
thors. Depending on the regularization parame- 
ter, this decrease is more or less pronounced (here 
we show the profile with the most pronounced de- 
crease). However, the ID inversion with the sec- 
toral splittings relies on the assumption that the 
rotation rate at any latitude is equal to the rate 
at the equator. 

As the differential rotation of the convective 
zone implies the higher latitudes rotate more 
slowly than the equator, we need to correct the 
splittings for this contribution. Therefore, we add 
12 nHz to the ^ = 1 splittings, 8 nHz to the £ = 2, 
and 6 nHz to the £ = 3, following Corbard et 
al. (1998). The values of these corrections were 
derived from a 2D Regularized Least-Squares in- 
version of MDI data in the convective zone, by 
comparing the splittings computed with the ex- 
act linear relation (that involves the 2D rotational 
kernels and rotation profile) with the splittings 
computed with the ID sectoral approximation. 
We present on the lower panel of Fig 2 the profile 
derived with the corrected splittings. We obtain 
a flat profile in the radiative interior. The dip 
around 0.3 Rq might be due to the absence of 
correction of the LOWL splittings but we cannot 
rule out a physical phenomenon at the limit of the 
nuclear core. 



This careful inversion of our most precise seis- 
mic data therefore favors a flat rotation curve 
at 430 nHz in the solar radiative interior. This 
lends support for a magnetic field strong enough 
to suppress any differential rotation that might 
arise from an angular momentum redistribution 
through the gravity waves (Talon, Kumar, & Zahn 
2002). We think that the absence of a decrease or 
increase in the rotation rate, compared with pre- 
vious works, is mainly due to the quality of the 
present data and the limitation in the degrada- 
tion of these data by the external layers. We note 
that an increase of the rotation rate seems now 
very unlikely down to 0.2 Rq, but cannot be ruled 
out below 0.2 by the present gravity-mode analysis 
(Turck-Chieze et al. 2002). 

We have also reduced the MDI-hGOLF split- 
tings to Oi coefficients (Schou et al. 1994) by re- 
moving 03 (for £ ^ 2,3) and 05 (for £ = 3) esti- 
mated from our knowledge of the rotation in the 
convection zone. This allows us to combine them 
directly with LOWL ai coefficients and carry a ID 
inversion for a latitudinally averaged rotation rate. 
This gives results nearly identical to the one shown 
on the bottom panel of Fig 2. Finally, we also sup- 
pressed the solar cycle effects for intermediate- and 
high-£ by using only the first year of LOWL data. 
Of course this increases the vertical error bars but 
does not change the rotation profile. 

6. Summary 

By using the high-quality seismic data from the 
spatial instruments GOLF & MDI, and by us- 
ing only the low-n modes, we limit the effects of 
the variable magnetic field that takes place in the 
outer layers and obtain a very coherent dataset for 
rotation below the solar convection zone. While a 
proper treatment of the latitudinal dependence of 
the rotation speed is still needed to extract the ro- 
tation in the core from ID inversions, the method 
we apply here removes the effect of latitudinal vari- 
ation of the rotation in the convection zone, and 
we obtain a fiat rotation profile down to 0.2 Rq. 
This puts a strong constraint on the redistribu- 
tion of the angular momentum. The uncertain- 
ties in the rotation rate are still quite large below 
0.3 i?o- progress on this point can be achieved by 
the detection of mixed low-£ pressure modes, and 
gravity modes. 
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Table 1 

Synodic sectoral splittings (Aa;„im) from GOLF-alone (left column) and MDI+GOLF 
(right column). The errors (err.) on the splittings are the la uncertainties provided by 
THE Hessian matrix of the fit. rt is the inner turning point of the mode. 
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Fig. 2. — Upper panel: rotation profile directly de- 
rived with GOLF & MDI, combined with LOWL 
splittings for £ > 3. Lower panel: rotation pro- 
file obtained after correction of the splittings to 
account for the differential rotation in the con- 
vective zone. The vertical error bars result from 
propagating splitting measurement errors through 
the inversion process, while the horizontal error 
bars give the FWHMs of the corresponding aver- 
aging kernels, and is an estimate of the resolution 
achieved at each depth. 
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